% This script generates Figure 2.5 of the notes
% Written by Dr Colin Caprani - www.colincaprani.com

xi = [0,0.05,0.1,0.5];
for i = 1:length(xi)
    [t u(i,:)] = sdof_damped(10,100,xi(i),0.025,0,4,0);
end
plot(t,u);
xlabel('Time (s)');
ylabel('Displacement (m)');
legend('Damping: 0%','Damping: 5%','Damping: 10%','Damping: 50%');